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A GENERAL MATHEMATICAL MODEL FOR BEAM AND PLATE VIBRATION 
IN BENDING MODES USING LUMPED PARAMETERS 

SUMMARY 


The analytic method proposed in this paper solves for natural frequencies and 
mode shapes of distributed mass and stiffness systems such as beams and plates. The 
structure is divided into an equivalent system of discrete masses and springs. A sche- 
matic diagram can be built up to show the nature of each item affecting the system. 
Increased off-hand evaluation of design changes results. A group of simultaneous 
equations is made up from the discrete schematic system. These equations comprise 
the mathematical model of the system. 

The mathematical model of the beam is completely general. The model is a 
matrix equation which will accept any input, any restraint, or any combination of end 
conditions. One model will suffice for all beam problems. Damping can be considered, 
but is not included in the main thesis of the paper. The model is in sufficient detail to 
accept the particulars of a nonuniform beam where mass and stiffness can be broken 
into lumped discrete but unequal values. The number of masses limits the eigenvalue 
solutions. Reasonable accuracy is dependent on a minimum of three masses between nodes. 
This is equivalent to four masses per modal number, plus one for the end node . Thus 
a system with seventeen masses would be reasonable for the calculation of the first four 
modes of a beam. 

Flat plates are simulated as two sets of cross-coupled beams. A nonuniform 
plate model is developed. End conditions of the individual beams may be varied, masses 
may be equated to zero, etc, to set up a rather complicated model. As in the beam model, 
the plate model is very general and only one model is necessary where the size is large 
enough to accept the required detail. Matrix size limits the number of eigenvalue 
solutions available. The model of a cross-coupled plate of sixteen masses is derived in 
the paper. Actually, practical solutions would require a much larger equation. 

The model does not account for any extensional deformation. Only bending is 
considered. The application of the method to shells and other curved plates is suggested 
when only bending modes need to be considered. 



SECTION I. INTRODUCTION 


Mechanical systems comprised of distributed mass and distributed stiffness have 
an infinite number of normal modes of vibration. Beams are characteristic of the distri- 
buted systems and have only one set of normal modes per single direction of deflection. 
They therefore are representative of the simplest approach to distributed system dynamic 
analysis. Other systems such as plates and shells can be represented by a system of 
cross-coupled beams once a suitable system of impedance and coordinate designation is 
selected. 

Systems comprised of uniformly distributed mass and stiffness lend themselves 
to the usual methods of solution, by utilizing some form of the Rayleigh method by equa- 
ting maximum strain to maximum kinetic energy for each mode . These methods arrive at 
mode shapes and natural frequencies but not finite response amplitudes. The incorpora- 
tion of damping into the calculations, then, allows an approach to finite amplitudes of 
response under steady-state conditions by equating input energy to dissipated energy. 

Some systems are not truly uniformly distributed but can be simulated by an equivalent 
uniformly distributed system. Typical of this class is a plate with equally spaced like 
stiffeners. The lower numbered modes of this type of stiffened plates will respond in 
close approximation to mode shapes of uniformly constructed plates. 

Some structures have irregular details of mass, stiffness, or end conditions that 
do not allow a convenient method of solution by distributed system analysis. They must 
be simulated by a system of lumped masses influenced by local spring effects. The use 
of these lumped methods is usually limited to a few of the lowest natural frequencies 
because a lumped system has a finite number of degrees of freedom and the deflection of 
each mass reflects a single point on the mode shape curve. Good definition of mode shape 
requires a minimum of about four masses between modes. Practical considerations of 
the problem solutions limits the number of useable lumped masses. Since design prob- 
lems are usually associated with low numbered modal response, the lumped-mass 
method could be of considerable importance. 

Any solution of a distributed system dynamic problem is of a complicated nature. 

A method that allows a simple subjective comprehension of each of the details has con- 
siderable value to the analyst. It will increase off-hand evaluation of the effect of any 
change, and lead to more effective design considerations. The analytic method proposed 
in this paper takes each specific detail and transforms it into an easily comprehendable 
part of the system diagram. The method is set up so that it is completely general. One 
program could be developed to accept all beam problems. Similarly, one program could 
handle all plate solutions. 

The method will be derived by detailed steps so that this paper can serve as a 
ready reference to users of dynamic analysis. 
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SECTION n. METHOD DERIVATION 


The first step in the method derivation is to review the basic equations relating 
moment to deflection (Fig. 1). When the deflection is very small and there is no net 



FIGURE 1. SECTION OF A BENT BEAM 

elongation of the beam, the deflection of the surface of the beam (BB', Fig. 1) can be 
considered parallel to the undeflected axis of the beam, and the deflection of the neutral 
axis of the beam (00’) may be considered perpendicular to the axis of the beam. If we 
arbitrarily set AL equal to OB, then simple geometric analysis will show that BB' equals 
00' in amplitude, and that they are 90 degrees apart in direction. The equation for the 
moment acting on section B'B' is 

^=Elff-. (1) 


where 


■M- = moment 
E = Young’s Modulus 
BB* 

A©= — -T— in radians for this case 
AL 

I = moment of inertia 
AL = unit length of beam. 

Now the moment is equal to the product of force, F, times arm length, AL. Equation 
(1) can be written 
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( 2 ) 


F AL = E I 


BB* 

(AL)- 


= E I 


00 ' 


(AL) 


or 

00 ' 

F = -^EI. (2a) 

Since the stiffness of the beam, K, is equal to F divided by the deflection 00', 

F El 

Equation (3) has the restrictions that the deflection be small and that moments acting on 
the stiffness must be acting through the arm length of AL. Figure I can now be simula- 
ted by a discriptive model (Fig. 2) . 



To simulate a beam. Figure 2 can be modified to represent a series of linked massless 
trusses with stiffness and mass lumped arbitrarily and with the basic concepts of Figure 
2 adhered to. Such a system is shown in Figure 3. It must be remembered that all 
effective motion of the spring attachment points is in the horizontal direction, and that 
of the hinges in the vertical direction (each related to the other). The interrelationship 




K = Spring constant 
M = Mass 
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FIGURE 3. SIMULATED BEAM 




of these two displacements leads to a defining of the spring deflection in terms of the 
deflection in the vertical axis. Each truss is hinged at the base to the adjacent trusses. 
As defined above, this hinge allows motion only in the vertical direction. Each hinge 
deflection is representative of the displacement of the neutral axis of the beam at its 
relative location. The relative displacement of one hinge of a truss to the other is equal 
in magnitude to the horizontal displacement of the spring attachment point. Thus all 
spring deflections can be equated to actual beam deflections. Let us define the net 


r‘ 

FIGURE 4. SPRING DEFLECTION MODEL 

deflection of spring A B, Figure 4, by the displacement of A by A^, etc. 

The following equation develops 

^A = ^2 - Ai, 

= ^3 ■ ^ 2 > 

^Ab"^b“^a"^3-2A2 + Ai. (4) 

The spring constants and the masses do not have to be equal for analytic solution. 
However the solution calculations are greatly simplified when all trusses are of equal 
size, except for the halves at the ends. This in effect divides the beam into segments 
of equal length between hinges, with the ends being half segments. Figure 5 depicts the 


(A) FREE END 


(B) HINGED END 


(C) FIXED END 
FIGURE 5. END CONDITIONS 
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three types of end conditions usually assumed in beam problems. The graphic descrip- 
tion is self evident. Now comes some of the real value of this method. Most deviations 
from the standard end conditions can be simply worked into the solution. Figure 6 shows 
a few of the possible non-standard end conditions of practical consideration. The non- 
standard impedance values attached to the ends are easily added to the general equations 
so that this type of problem is included in the general solution. 



^77X77777777777777^ 



(D) RESTRAINED HINGE 


(E) FLEXIBLE ATTACHMENT 


(F) END MASS RESTRAINT 


FIGURE 6. NON-STANDARD END CONDITIONS 


Attachments to the beam also can be accounted for in the general solution. Fig- 
ure 7 shows the diagramatic inclusion of such attachments to the simulated system. 



(G) FIXED ATTACHMENT 




(H) SPRING ATTACHMENT (I) MASS ATTACHMENT 


FIGURE 7. BEAM ATTACHMENTS 
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It should be noted that damping has not yet been inserted into the diagram. A 
type of damping simulation is possible but not very practical within the present state-of- 
the-art. Its usefulness will depend upon further research. The diagramatic model de- 
scribed up to this point is for free vibration systems. The characteristics of the system 
can be conveniently represented by impedance values. Such representation can lead to 
an understanding of the damping properties of the system. Structural damping is inserted 
into the impedance equation as follows: 

Z = jwM + 


• IVT ^ 

icoM + ' . 

JW O) 


( 5 ) 


M = mass 
K = stiffness 
rj = damping loss factor 
(jo = angular frequency. 

Force, F, is related to displacement, A, by the relation 

F = (Z) (jwA). (6) 

The simulated system shown in Figure 3 can be solved by using the usual equations of con- 
tinuity of force and moment. Damping can be included in the impedance equation. How- 
ever when used, it is recommended to solve the reactive portion first without damping, 
and then insert damping to calculate amplitudes. For simple beam and plate models, 
damping can be schematically added as in Figure 8. Any of the usual damping concepts 
are useable but here we use a damping symbol that is indicative of slipping or Coulombs* 



damping and this type is proportionate to displacement. Thus the damping is proportionate 
to the same displacement that works on the spring located between each of the trusses. 
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The loss factor is unit length or area dependent. Once the mode shape is determined 
for a mode the input energy must equal the output energy, 



= EftKAj 


N 


(7) 


where m is representative of forced input locations and n is representative of response 
deflections of the various springs. Equation (7) would be applicable to beams or plates 
without attachments. By including representative attachment loss factor increments, 
equation (7) could become 



E 

N 


jT? K A2 + 


2 K a2 

N» 


( 8 ) 


where N*, 77 ' and A^^, refer to the damping caused by the attachments. The schematic 
inclusion of the damping caused by the attachments is represented in Figure 9. Equation 
(8) will solve for discrete values of amplitude for each mode. The A deflections are 
dependent on the y displacements. This aspect will be discussed after development of 
the general solution equation. 



FIGURE 9. INCLUSION OF ATTACHMENT DAMPING 

It has been the purpose here to show the possibility of damping analysis by using 
the analytical model. The remainder of the report will deal with undamped systems only. 


SECTION m. BEAM SOLUTIONS 


Development of a General Beam Matrix Equation 

The beam solution uses a set of equations that define the continuity of force and of 
moments. Figure 10 is a typical model of a beam of six elements of mass. All motion 
of the beam is assumed in the ±y directions only. End restraints may impose a reaction 
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FIGURE 10. MODEL OF A SIX MASS BEAM 


FREE 

■Mo = 0, Ro = 0 

FIXED 

yo = 0. yi = 0 

HINGED 

yo = 0, M^j = 0 


FREE 

~ 0, Ry — 0 

FIXED 

ye = 0. y? = 0 

HINGED 

Y7 = 0, -My = 0 


Mo, M 7 are end moments, in. lbs. 

Ro, R 7 are end support forces, lbs. 

El 

Kj etc. is stiffness constant -Jz> lbs. /in. 

Ai etc. is the net deflection of the spring "K”, in. 
yo etc. is the neutral axis deflection of the beam, in. 
i is the height and base of each full triangle , in. 

Ml etc, is the lumped mass of the section of the beam. 


force R and/or a moment ofM; The end may have a displacement of yo or y-j. Each 
element of mass M and stiffness K may have its own assigned quantitative value. In the 
case of a uniformly distributed system, all M values would be equal and all K values 
would be equal. 


The Symbols used in the derivation of the equations are defined below Figure 10. 
The counter-clockwise momentM-o is considered positive, whereas the clockwise mo- 
ment M 7 is considered positive. This is the usual convention used in beam equations. 
The deflection of y in the downward direction is considered positive and the deflections 
of the springs to the right are considered positive. Both support forces, Rq and R 7 , are 
considered positive in the upward direction. 
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The first set of equations to be evolved are based on the reasoning leading to equa- 
tion (4) . We define Aj as the deflection imposed on and A 2 as the deflection imposed 
upon K 2 , etc. , a positive value of A reflects an elongation of the spring. Bemember 
that the moment arm of the springs is equal in length to the distance between hinges. 


-Ai + y2 - 3yi + 2yo 

(9) 

-A 2 + ya - 2y2 + yi 

(10) 

--^^3 + 74 - 2ya + ya 

(11) 

-A 4 + yg - 2y4 + ya 

(12) 

-Ag + yg - 2yg + y 4 

(13) 

-Ag + 2y7 - Syg + yg 

(14) 


In equations 9 through 14 the y’s are actual deflections at the locations as indicated by 
the subscript. 

Next, in Figure 11, let us look at the moments relative to a section 1, 1 of the 
beam passing through location 1. The deflection in spring is induced by equal and 
opposite moments on either side of the section 1, 1. For convenience let us take the 
moments to the left. Then the continuity of moments at section 1, 1 will be: 


Aj Ki f - •M’q - 2 Rfl- 

This can be resolved to a more standard form of: 

^0 

0 — Aj Kj - ^ + . 5 Rq. ( 15) 

Note that-Mo and Rq result in moments of opposite sign as tiiey resist each other. 

Section 2,2, Figure 12, now involves more than end conditions in the moment 
equation. Since the springs at locations other than the designated section location only 
transmit moment, they are not involved in the equation, "hto and Ro, as external effects, 
do not represent physical characteristics of the beam. The mass, on the other hand, 
is an impedance value and can be treated as such. The moment about section 2,2 in- 
duced by any impedance, Zj, with a net displacement of y^ would be: 
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FIGURE 12. BEAM SECTION AT LOCATION 2 
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^(Zl) = (Zi) (j wyj) i. 


For Zj in this case we can substitute j co then: -M(Mi)= - Mj yi i . The moments 
at section 2,2 can now be equated as: 

iVt 

0 = Ag Kg - + 1. 5 Ro - 0)2 yi Ml. (16) 

Following the same reasoning as above, equations of moment continuity can be written 
for sections at 3, 4, 5, etc. For a six mass beam the set of moment equations would be: 


^0 

0 = Ai Ki - -j®- + . 5 Ro 




0 = A 2 K 2 - + f. 5 Ro - 0)2 yi Ml 


(15) 

(16) 




0 = Ag Kg - "Y" + 2. 5 Ro - 20)2 Ml - o)=^ y 2 M 2 

0 = A 4 K 4 - -p + 3 . 5 Ro - 30)2 yi Ml - 2o)2 yg M 2 - 0)2 yg Mg 


(17) 

(18) 


M- 

0 = Ag Kg - + 4. 5 Ro - 40)2 y^ Mi - 30)2 Mg - 2o)2 yg Mg - 0)2 y 4 M 4 


(19) 


0 = Ag Ke - + 5. 5 Ro - 50)2 yi Mi - 40)2 yg Mg - 30)2 yg Mg - 2o)2 y 4 M 4 

- 0)2 yg Mg 


( 20 ) 


We can also write the moment equation at y^ and it becomes: 


0 



+ 6 Ro - 5. 5 0)2 yi Ml - 4. 5 0)2 yg Mg - 3. 5 0)2 yg Mg 


- 2 . 5 0)2 y^ M 4 - 1. 5 0)2 yg Mg - 0. 5 0)2 yg Mg. 


( 21 ) 


If we substitute equation 9 through 14 into equations 15 through 21 as applicable, 
we will end with the following list of unknowns: 

yo, ^0. Ro» yi. Y2> ys» y4. Y5> ye» y?.^* R?- 
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Any choice of end conditions will reduce this group by four, which gives eight unknowns. 
Thus far we have developed seven equations, numbered 15 through 21. We need one 
more to satisfy our simultaneous equation solution. In forced vibration the sum of the 
forces acting in the vertical direction must equal zero. The vertical forces are induced 
by the end forces Rq and R 7 , and the inertial forces by the masses or any other restraint 
at the hinges. The eighth equation can now be developed: 


0 = - Ro - R7 + yi Ml + 0)2 yj M2 + 0)2 yg M 3 + 0)2 M 4 + 0)2 yg M 5 + 0)2 yg Mg 

( 22 ) 


Combining equations 15 through 22 produces matrix equation (23). 


0 = 


.5 

0 -i'* 

0 

+2^1 

1 

CO 

-Ki 

0 

0 

0 

0 

1.5 

-.1 

0 

0 

0 

+ (K2-W%i) 

-2K2 

+K2 

0 

0 

0 

2.5 

0 

0 

0 


+ (K 3 -o; 2 M 2 ) 

- 2 Kg 

+K3 

0 

0 

3.5 

0 

0 

0 

- 3 u%i 


+(K 4 -w 2 m 3 ) 

-2K4 

+K4 

0 

4.5 

0 

0 

0 



-2oj^M3 

+ (K 5 -w 2 M 4 ) 

- 2 Kg 


5.5 

0 

0 

0 




-2U)2m4 

+ (Ke- 0 ) 2 M 5 ) 

- 3 Ke 

6 

0 

+i' 

0 

- 5 . 

- 4 . 

- 3 . SW^Mg 

-2. 5 w^M 4 

-1. 

-. 50)2 

-1 

-1 0 

0 

0 



+oi^M3 

+w 2 M 4 

+ 0)%5 

+w2M( 


Matrix Equation, Beam General Model of Six Masses 


Equation (23) 


Ho| 

R 7 

Mo: 

M 7 

yo 

yi 

yz 

yz 

y4 

y% 

y% 

y? I 


Equation (23) represents a completely general matnematical model of a beam. 
The rectangular matrix and the right-hand column matrix will reduce to solvable square 
and column matrices of the same order upon selection of the end conditions. Note, from 
the table below Figure 10, that each end condition has two of the possible variables equal 
to zero. The selection of end conditions would then reduce the rectangular matrix of 
equation (23) to a square 8 x 8 matrix and the right-hand column matrix will have four 
of its variables equal to zero and therefore cancel out to the 8 th order. 
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Since any selection of end conditions will reduce the variables by two , the end 
conditions need not be alike. The matrix equation will accept any combination of end 
conditions. This is a distinct advantage of this type of approach. Once the computer is 
programed to solve the Eigenvalue and Eigenvector problems of full, non-symmetrical, 
matrices the analytical model will produce solutions not convenient by other methods. 
Since the matrix is not dependent on uniform constants, a non-symmetrical physical con- 
figuration may be solved. In a formal sense of application this method would apply only 
to stepped structures with uniform constants between steps. In a loose sense the triang- 
ular segments may represent a lumping of the variable constants into equivalent values. 
For instance if the mass of the beam is equatable to a function of position on the beam. 


M 

n 


X 


(X) 


n-1 


dx 


where = Mass of the nth segment of the beam 

p = Mass density of the beam in mass per unit length 
X = distance along the neutral axis of the beam. 

The deflection of the beam must be considered linear and estimates made must be 
for the spring constant of each segment of an irregular beam. 

Equation (24)shows the reduction of equation (23) to represent a uniform symme- 
trical beam of six masses. This specific model is used to illustrate the 3 usual beam 
conditions of 


1. Fixed-Fixed 

2. Hinged-Hinged 

3. Free-Free 


This simplified model, equation (24) , was used because it 
-£ +2K 


0 = 


.5 

1.5 
2. 5 

3.5 


-1 


-1 


-i 


-3K 

(K - w^M) 


+K 
-2K 
+(K - 

+ (K - 2oj^M) 


+K 

-K 

■(K + w^M) 


1^0,7 
J*^,7 
yo,7 
yi 6 


Matrix Equation, Symmetrical Uniform Beam General Model of 6 
Masses Equation (24) 


ya.s 

ys,4 
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allowed reasonable (4 by 4) hand solution to show examples of problem solution. Higher 
order solutions are recommended for the practical use of the mathematical model. 

The solution of matrix equations such as equation (23) is somewhat involved. 
When the rectangular matrix is reduced to a square matrix it represents the impedance 
of the system. At each resonant frequency the reactive portion of the impedance is equal 
to zero. Since only reactive impedance has been inserted into the rectangular matrix, 
the values of co that equate it to zero are solutions of the natural frequencies. These 
values of co are known as Eigenvalues. 

The variables in the column matrix represent system response. Each natural 
frequency will have a finite set of response values. Solutions of the column matrix are 
known as Eigenvectors. Even after the insertion of a natural frequency into the square 
matrix, equation (23) is not solvable. However, if we normalize the response values to 
one of the Eigenvectors, e.g. , yi, then the insertion of unity into the column 
matrix makes the equation solvable. This means that only relative displacements can 
be determined from an Eigenvalue - Eigenvector solution. These relative values de- 
scribe a mode shape for each natural frequency. 

The physical reason finite values of response amplitude cannot be had from a 
solution of equation (25) is that response amplitudes of a steady state forced system 
are damping dependent. Equation (25) contains no damping values. 

Equation (8) equates input energy to dissipated energy. Each deflection y is 
equated as 


y = K y , 
n n-^y 


(25) 


Where y^, is the normalization factor. All K^'s must be determined. Then the substitu- 
tion of equation (25) into equation (9) through (14) and the resultants into equation (8) will 
lead to solutions of finite deflection values for each mode. Usually, at the first few of the 
lowest natural frequencies such as those suitable for this method of analysis, these de- 
flections will be representative of total deflections at the natural frequencies. This is 
true because modes other than the resonant mode will have very little amplitude when 
only low frequencies are being considered. 


Examples -4x4 Matrix Eigenvalue Solutions of Uniform Beams 

The general solution, equation (23), can very often be reduced in scope. This 
is especially true when the beam is symmetrical and uniform. In fact, equation (24) is 
the first step in reduction to account for uniformity. When the beam is symmetrical it 
has an even number of masses and the odd modes only are considered. Each deflection 
on the right of center should have an equal deflection on the left side. Also the moment 
and reaction force at the ends are equal. By substitution of one set into the other. 
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equation (24) can be reduced to half size. Equations 26, 27, and 28 represent uniform 
symmetrical beams of six masses that reduce to a 4 by 4 system. These 4 by 4 matrices 
are easily calculated by hand and will be used to check out the accuracy of this analytic 
method ly general comparison with standard solutions. 


Matrix Equations, Symmetrical Uniform Beam, Specific CaSes 


1. Fixed - Fixed: 


.5 


+K 

0 


Ro^7 

1.5 


-2K 

+K 


Moj 

2. 5 


+ (K - w^M) 

-K 


5 

3. 5 


+(K - 2o)2m) 

-(K + w^M) 


ya 4 


2. Hinged - Hinged : 


.5 

-3K 

+K 

0 


Rq^ 7 

1. 5 

+(K - w^M) 

-2K 

+K 


yi,6 

2. 5 

-2(Jm 

+(K - co^M) 

-K 


y2,5 

3.5 

-3co2m 

+(K - 205%) 

-(K + 0)%) 


ys.4 


3. Free - Free: 


2K 

-3K 

+K 

0 


yo,7 

0 

(K - CD^M) 

-2K 

+K 


yi,6 

0 

-2o)2m 

+(K - w^M) 

-K 


y2,5 

0 

-3w^M 

+(K - 2<Jyi) 

-(K + OJ^M) 


y3,4 


(26) 


(27) 


( 28 ) 
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Equation (29) is a longhand solution of a four by four matrix. The calculation 
will be based upon this solution in tabular form. It should be noted that when A, B, C 
or D is equal to zero, a portion of the tabulation need not be performed. For example, 
D equals zero in the solution of the fixed-fixed beam, equation (25). 


4 by 4 Matrix Expansion Model 

A B C D 
E F G H 
I J K L 
M N P Q 

= AFKQ + AGLN + AHJP - AFPL - AJGQ - ANKH 

- BEKQ - BGLM - BHIP + BEPL + BIGQ + BMKH 
+ CEJQ + CFLM + CHIN - CENL - CIFQ - CMJH 

- DE JP - DFKM - DGIN + DENK + DIFP + DMJG 


TABLE I. MATBIX EXPANSION OF EQUATION (26) 



K^ 

Kw^m oj^M^ 

AFKQ = . 5 (-f “^) (K - OJ^M) (-K - w^M) 

= + . 5 f K^ 

-. 5f 

AGLN = . 5 ( -2K) ( -K) ( -i “^) 

= K^ 


AHJP = . 5 (K) (-i”^ (K - 2 o;2m) 

= -.5 f"* K^ 


-AFPL= 5(-f "^) (K - 2^2 m) (-K) 

= -.5f K^ 

+f w^mK 

-AJGQ = 5( -1 ( -2K) ( -K -w^M) 

= +f k2 

+f 0)2MK 

-ANKH = 5(-i"b (K - w^M) (K) 

= + . 5i K^ 

- . 5 f w^MK 

-BEKQ=+f~^ (1. 5)(K - cj2m)(-K - w^M) 

= - 1.5f K^ 

+ 1. 

-BGLM = + f “^ ( -2K) ( -K) ( 3. 5) 

= 4- 7 f K^ 



= A 


F G H 


E G H 


E F H 


E F G 

J K L 
N P Q 

- B 

I K L 
MP Q 

+ C 

I J L 
MN Q 

- D 

I J K 
MN P 


(29) 
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Table I (Cont'd) 


-BHIP = + 1 (K) (2. 5) (K - 

= + 2. 5i K2 

- 5 i 

+BEPL= -i"* (1.5)(K - 2 o)2m)(-K) 

= + 1.5i"^K2 

- 3f“^MK 

+BIGQ = (2. 5) (-2K) (-K 

= -5i 

- S 

+BMKH = - i ”^3. 5) (K - (JU) (K) 

= - 3. 5 « 

+ 3.5i"^a;2MK 

+CEJQ = +K( 1. 5) ( -SL ■*) ( -K - co^M) 

= + 1. 5i~^2 

+ 1.5i"VMK 

+CFLM = +K(-f "^) (-K) (3. 5) 

= + 3. 5i 


+CHIN = +K(K) (2.5) (-i"b 

= - 2. 


-CENL= -K(1.5)(-«"^)(-K) 

= - l.Si'V 


-CIFQ = -K(2.5) (-je"^ (-K - co^M) 

= - 2.5f 

- 2. 5 f Cc)2]y[j^ 

-CMJH= -K(3. 5) (-i (K) 

= + 3. 5i 



4. Fixed-Fixed Beam Sample Solution 


Equation (26) represents the solution of a uniform symmetrical fixed-fixed ended 
beam. The square matrix equals zero for any natural frequency that can be calculated 
from the mathematical model. First the matrix has to be expanded to a single equation. 
By using equation (29) as a model, Table I is constructed. The summation of Table I 
set to zero is 

0^ - 8^2 M K + SK^ = 0. 

If we use the quadratic solution equation of 


18 





9 -b ± 
(jO^ = 


4^ 


4 ac 


2a 


and set 


a = M® 
b = - 8 M K 
c = 3 K^, 

then CO* can be solved as follows: 


co* = 


+ 8 M K ± 64 M* K* - 12 M* K* 


2 M* 


9 + 8 M K ±'s/s2 M* K* 


co‘ = 


2M“ 


9 + 8 M K ± 7. 2111 M K 

= 5 

2 M* 


9 . 394 K 7. 605 K 

co‘ = T 

M ’ M 


Since we have shown K 


El 




CO 


2 

1 


. 394 E I 
M i* 


The standard equation for the first natural frequency of a fixed-fixed beam is 


CO 


2 

1 


(22.4)* E I 

Mt L* 


In this equation M^. = 6 M and L = 6i. These, when substituted in, yield 



386 E I 
M i* 
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The error in using the mathematical model of equation (26) would be 


Error = 


394 - 386 
386 


X 100 = 2. 1%, 


when o)^ is considered. Actually the error in terms of u> would be 


100 X 


's/3M -\l386 
n/386 


= 100 X 


19. 85 - 19. 65 
19. 65 


20 

19. 65 


1. 04%. 


Since the second value of as calculated is about 20 times the first, it is obvious that 
it does not represent the third mode. We can conclude that the four by four matrix can 
be good only for the one lowest natural frequency calculation. If an eight by eight ma- 
trix were used the accuracy of the first mode would be increased and the third mode 
made available. Figure 13 is the diagram for equation (26). It is noticeable that only 
two masses are in oscillation. This is not sufficient freedom to simulate a good third 



FIGURE 13. FOLDED MODEL FOR FIXED-FIXED BEAM 

mode. That is why the third mode can not be calculated by this small size matrix. 

It should be remembered that the symmetrical method eliminates the second mode. 

For the first mode the column matrix of equation ( 26) can not be determined in 
absolute values. A solution can be performed in the following manner by using equation 
(26): 

Substitute in the calculated value of 
Divide the column matrix by y^. 
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Solve the set of simultaneous equations for R/ys, ■M(/y 3 and y^/Yz- 

Note that R(/y 3 is the reaction force for the first mode at the beam attachment location 
per unit deflection of the beam at y^ and that-M(/y 3 is the moment resulting from the 
first mode at the beam attachment per unit deflection at y 3 . 

Other aspects of the system can be calculated. The average slope of the beam 
for the first mode between y 2 and y 3 is 



in radians per unit deflection at y 3 . The maximum stress in the beam is 

S = — , (30) 


where C is the distance from the neutral axis to the extreme fiber. Since we have cal- 
culated -Wys we get 


_S 

ys 


y3 1 


(31) 


as the maximum stress per unit deflection of ys. It must be remembered that stress 
calculated by equation (31) is for one mode of vibration. The total normalized stress 
resulting from all modes would be the complex summation of 



where Sj^/y Is the normalized maximum stress of the ntn mode. 

Hi nged-Hinged Beam Sample Solution 

Equation 27 is a matrix solution for a uniform symmetrical beam with both ends 
hinged, or simply supported. Equation 27 is solved using equation 29 as a model. Table 
n is a tabulation of the calculations. 
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table II, EQUATION (21) MATRIX EXPANSION 


AFKQ = . 5(K - w^M) (K - co^M) (-K - w^M) 

,5o;®m3 +. 5o;%% +. 5K® 

AGLN = . 5(-2K) (-K) (-30;%) 


-3.0 


AHJP = . 5(K) (-2co%) (K -2co%) 

+2. 0 

-1.0 


-AFPL = 5(K - w^M) (K - 2co^M) (-K) 

+1.0 

-1. 5 

+. 5 

-AJGQ = 5(-2cu%) (-2K) (-K - cu%) 

+2. 0 

+2. 0 


-ANKH = 5(-3o;%) (K - co%) (K) 

-1.5 

+1. 5 


-BEKQ = 3K( 1. 5) (K - cu%) (-K - a;%) 

4.5 


-4.5 

-BGLM = 3K(-2K) (-K) (3. 5) 



+21.0 

-BHIP = 3K(K) (2. 5( (K - 2o)%) 


-15. 0 

+7. 5 

+BEPL = -3K( 1. 5) (K - 2w%) ( -K) 


-9. 0 

+4. 5 

+BIGQ = -3K(2. 5) (-2K) (-K - o)%) 


-15. 

-15.0 

+BMKH = -3K(3. 5) (K - co%) (K) 


+10. 5 

-10.5 

+CEJQ = K(1.5) (-2cu%) (-K - o)%) 

+3.0 

+3, 0 


+CFLM = K(K - u;%) (-K) (3. 5) 


+3. 5 

-3. 5 

+CfflN = K(K) (2. 5) (-3co%) 


-7.5 


-CENL = -K( 1. 5) (-3c<;%) (-K) 


-4.5 


-CIFQ = -K(2. 5) (K - 0 ;%) (-K - u;%) 

-2.5 


+2. 5 

-CMJH = -K( 3. 5) ( -2o;%) (K) 


+7. 0 



The summation of Table II yields: 

0 = . _ 28. Sw^MK^ + 2K^ 
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rr i 2 X K 

If we set = — — 
M 


the above equation becomes 


0 = . 5 X® + 9 - 28. 5 X + 2. 

We know that x will be a small number in the order of 0. 1. Therefore the x^ term in the 
equation can be neglected. Use the quadratic equation general solution to get 

18 

then 

> , 0718 K . 0718 E I 

00 ^ = = 5 

M M ■ 

The equation for a uniform symmetrical hinged beam is 


CiJ 


2 

1 


. 0752 E I 
M 


The error for is 


100 X 


. 0752 - . 0718 
. 0752 


The error for co is 


4. 5%. 


100 X 


'sp752 - 's/718 
nT752 


2. 3%. 


As stated in the fixed-fixed beam example, the 4 by 4 matrix is in reality too 
small for good problem solution. The third mode solution is not available. The normal- 
ized moment, end resultant force and deflections for the first mode can be calculated as 
well as the average slope between hinges. 

TABLE III, MATRIX EXPANSION OF EQUATION ( 28) 


2 K(K-o;2 M) (K-co^ (-K-w^M) 

2 K(-2K) (-K) (-30)2 jyjj 
2 K(-2 o)2 M) (K -2o)2 M) (+K) 


-2K^ + 2K2o)2 M + 2R2 oj^ - 2Ko)®m2 

- 12 K® 0)2 M 

- 4K2 0)2 M + 8 k2 o)^ m2 
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-2 K(-3w^ M) (K - 0)2 M) (K) 

-2 K(K-o)2 M) (K -2o)2m) (-K) 

-2 K("2 0)2 M) (-2 K) (-K -0)2 M) 


Table IH (Cont'd) 

= 6 k2 0)2 M - 6 k 2 o)^ M* 

= 2 - 6 k2 0)2 M + 4 k2 0)^ m2 

= + 8 0)2 M + 8 0)^ m2 


- 6 K® 0)2 M + 162 k 2 o)^ m2 - 2 K O)® M® 


Free-Free Beam Sa mple Solution 


Equation ( 28) is the matrix solution for a uniform symmetrical beam with both 
ends free. The square matrix of equation (28) reduces to: 


2 K 


K - 0)2 M 

- 2o)2 m 

- 3 0)2 M 


- 2 K 

+ K - 0)2 M 
+ K -2o)2 M 


K 

K 

K - 0)2 M 


The solution of this matrix when set equal to zero is from Table HI and is 
0 = - 2 o)® M® K + 16 0 )^ m2 k2 - 6 0)2 M K®. 

X K 

If we let 0)2 = — ^ , then the above equation reduces to 
0 = + 2 x2 - 16 x+6 = . 395. 

Then, 

9 _ . 395 K .395 E I 
‘*^1 M i®M ■ 

The first natural frequency of the free-free beam should be 


OJ 


2 = 
1 


. 386 E I 
Mi® 


The error in terms of o)® will be 
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100 X 


. 395 - .386 
.395 


2. 3%. 


The error in terms of co will be 


100 X 


V395V386 

'v/396 


1 . 1 %. 


Ca ntilever Beam Sample Solution 

The cantilever beam is a special case of the beam equation. It has one end fixed 
and one end free. Thus it is not a symmetrical beam. Figure 14 represents a diagra- 
matic model that fits a 4 by 4 matrix. This particular model is used because it utilizes 
the same model configurations of equations 26 through 28. Moments are taken at the 
hinge locations. Equation (3 2) represents the summation of moments about the first 
mass, M^, by using the left input to balance the moment of spring Kj. Equation (33) lists 
the summation of moments that are to the left of mass M 2 and reacting on spring K 2 . 
Equation (34) represents the summation of the moments acting to the right of mass M 2 
and acting on spring K 2 . This switching of sides is convenient because spring Kj, ac- 
cording to this model, has no deflection. Equation (35) is the summation of forces acting 
in the vertical direction. Equations 36 and 37 equate the spring deflections to hinge dis- 
placements. Equation 38 follows from equations (32) through (38). Table IV is the 
solution of the square matrix of equation (38). 


A j A2 A3 



FIGURE 14. MODEL OF A CANTILEVER BEAM 


25 



Cantilever equations 


0 = -AiKi+-^ - .5Ro 

0 = -A2K2 + ^- 1.5Ro 

0 = -A2K2 + 


0 = -Ro + w^M2y2 


=Y2 

^2 = ys - 2yj 


i-1 

-.5 

-K 

0 


itf-0 

i"^ 

-1.5 

+2K 

-K 


Ro 

0 

0 

+2K 

-K + co*M 


Y2 

0 

-1 

co*M 

co*M 


Y3 


(32) 

(33) 

(34) 

(35) 

(36) 

(37) 


(38) 


TABLE IV. MATRIX EXPANSION OF EQUATION (38) 


(-1.5) (2K) (Co^M) 

(2K) (-1) (-K + w^M) 
i‘l (0) 

(-1) (2K) (-K) 

-i"^ (-1.5) (o)2m) (-K + co^M) 
-i-i (0) 

.5 (i*^) (+2 K) (UJ^M) 

.5 (0) 


= -3co^KMi-^ 

= 2K^ -2u)2KMi"^ 

= 0 

= -2 

= -1.5 CO* K M 

= 0 

= +co*KMi"^ 

= 0 


+1.5 co^ M*i-^ 
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Table IV (Cont'd) 


.5 (0) 

-.5 (0) 

-.5 (0) 

5 £~^(co^ M) (-K + 0)2 M) 
-K (0) 

-K (0) 

-K (0) 

+K (0) 

K (0) 

K i"^ (-1) (-K + 0)2 M) 


0 

0 

0 

+0.5o)2KMi"^ -0. 5o)^M2i-^ 

0 
0 
0 
0 
0 

K2i-1 -0)2 KMi"^ 


Total of Table IV 

Setting the square matrix equal to zero we get 
0 = O)^ m 2 - 6 o;2 K M + k2. 


Let 


O) 


2 


X K 
M 


then 


0 = x2-6x+1, 


X = 


6 ±'s/36 - 4 
2 


6 ± 4 ^ 
2 


-6 0)2 K M +o)^M2i“^ 


X = 


6-5.67 

2 


.33 

2 


then 


. 165, 
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CO 

1 


. 165 K _ . 165 E I 
M j?® M * 


The standard solution for the first mode of a cantilever beam will give 


CO 


2 

1 


. 153 E I 

M ■ 


The error in terms of co^ is 


. 165 - . 153 
. 153 


X 100 = 7. 9% . 


The error in terms of co is 


Vl53 


X 100 = 3. 2%. 


Since there are only two displacement values in this example, and ys, good accuracy 
can not be expected. 


Summary of the Worked Examples 


The four examples worked here are algebraic with the end resultant in a form 
that could be compared to the general beam natural frequency equation. For an example 
where the masses, stiffnesses, and length dimensions can be determined, the Eigenvalue 
solution of the square matrix will find the natural frequencies directly without the substi- 
E I 

tution of -jy- in for K at the end of the problem. Actually this substitution will be made 


when determining the initial values of stiffness to insert into the matrix. It need only be 
made once. 


Although a 4 by 4 matrix is obviously too small to assure enough Input detail, 
the accuracy of the four problems is good enough to point out the validity of the method. 


Non-Uniform Beam Solutions 

The mathematical m'odel put forth in this paper is adaptable to beams of non- 
uniform stiffness and/or mass as well as beams with a curved neutral axis. These prob- 
lems can be developed to fit a standard Eigenvalue solution to a square matrix but the 


28 


matrix itself is, of course, not standard in form. Each of these special situations must 
have a special set of equations developed. First we must examine the limitations of the 
method. For an example, consider Figure 15 as typical of a special beam base. 



FIGURE 15. NON-UNIFORM, NON SYMMETRICAL BEAM 

Each link of the simulated model. Figure 16, must have like dimensions but the indivi- 
dual values of M and K can vary. The trussed links must be small enough so that the 
net deflection of the beam hinge can be considered a net deflection of the spring. Error 
will be proportionate to the relative angle mode by the bases of adjacent trusses. This 
angle will increase as the length "i" is increased, and decrease as the length "i" is de- 
creased. The allowable error puts a limitation on the model truss size and therefore 
on the minimum number of trusses. 



FIGURE 16. NON-UNIFORM BEAM MODEL 

The deflection at M^ will have a resultant with both an "x” and a "y” component. 
However, if a single coordinate system of response relative to location is worked out, 
a smaller matrix results. 
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In this particular model there is no moment on spring K^. The motion of each 
mass concentrated at the hinge is considered to be in the direction of a ray that bisects 
the angle between the adjacent truss bases. The moment arm about Mj will then be a 
perpendicular struck through the location of to the directional ray. By using this 
method a set of moment equations can be developed in the same manner as for the straight 
beam. The force summation must equal zero in both the x and the y direction. A set of 
equations can be developed to set up a matrix equation that will solve for Eigenvalues and 
Eigenvectors. 

The use of the mathematical model presented herein probably has its greatest 
value in helping to solve the more complicated situations by reducing the specifics to 
a model that can be solved. 

SECTION IV. VARIATIONS ON THE ANALYTICAL BEAM MODEL 


Initial Displacement 

If an initial displacement is known, it replaces the appropriate y. For example: 
if y 4 has an initial displacement of .005 inches, this value will replace y 4 in the column 
matrix of equation (23) . 


Mode Shape 

The maximum displacement is given an initial displacement of unity. A plot 
of the values of y versus location on the beam will describe the mode shape. Curve 
fitting procedure will equate y to x. 

Forced Oscillation 

Forced oscillations will enter the matrix equation (23) as a column matrix on the 
left side. The procedure is illustrated by Figure 17. 



Fi 


F2 


FIGURE 17. FORCED EXCITATION EXAMPLE 
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0 


M- 

= - AiKi+-p-. 5 R^ 

+ ,5Fi A 2 K 2 + ^ - 1. 5 Rg 

M- 

+ I.5F1 = - A3K3 - 2.5 R(, + 0)2 MgYg 

M- 

+ 2. 5 Fi = - A 4 K 4 + -3.5 Ro + 2o)2 Mj M ^3 

+ . 5 F 2 ~ “ ^4^4 

+ (Fj + F2) = - Ro + M2y2 + ^2 M^2 ^373 + M4Y4 

Ai = Y2 
A 2 = Y3 - 2 Y 2 
A3 = Y4 - 2yg + 

A 4 = 2 y^ - 3y^ + yg 


0 


- . 

5 

+i-l 

-Ki 

0 

0 

0 


R 

0 

+ . 5 F 4 


-1. 

5 

+f-l 

+2K2 

-K2 

0 

0 



+ 1. 5 Fi 


-2. 

5 

+i-l 

(-K3 + 0)2 M2) 

+2K3 

-K3 

0 


^2 

+ 2. 5 Fi 


- 3 . 

5 

+i"^ 

+ 2 o )2 M 

(-K4 + 0)2 M) 

+ 3K4 

-2K4 


^3 

+ .5F2 



0 

0 

0 

-K4 

+3K4 

-2K4 


^4 

+ (F4 + F2) 


- 

1 

0 

0)2M2 

o)^M3 

o) 2 M 4 

0 


^5 


Mixed End Conditions 

When the end of the beam is anchored to a flexible member, f A© Kq is sub- 
stituted in equation ( 23 ) for Mo and Aq — 2 y^; y^ is then a variable and replaces ^ 4 o ^ 
right side column matrix. 
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When the end condition has a rotational oscillation, fCx^) is then a known input 
f(x ) = 2yi and its function is substituted into the right hand column matrix for yj. 

When the end is oscillating in the y direction, f (y©) is a known function and is 
substituted into the right hand column matrix (Fig. 18) . 



Multi-Supported Beams 

For a multi supported beam, must eliminate y^ in the matrix equation. Thus 
Rjj will replace y^ in the basic equations and R^^ will replace yj^ in the right hand 

column matrix as an unlmown (Fig. 19) . 



FIGURE 19. EXAMPLE OF MULU SUPPORTED BEAM 
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In the non-rigid supported beam, as illustrated, the support must be at a hinge 
between triangles. In each instance of support, for example where and 

Mg yg are in the matrix (23) , the quantities (co^ - K^) and Mg - Kg) are re- 
spectively substituted since the spring will resist motion at the support. If the springs 
are initially not at rest, the respective forces must be inserted (Fig. 20) . 


K K K K K K K 



FIGURE 20. EXAMPLE OF FLEXIBLE SUPPORTED BEAM 

SECTION V. ADAPTION OF A BEAM ANALYTICAL MODEL 
TO A PLATE ANALYTICAL MODEL 


The beam of Figure 3 is simulated by a model that is most convenient for all of 
the beam problems. Other models may be constructed to meet specific purposes. The 
beam of Figure 3 is not readily adaptable to plate configurations. Therefore, it is 
necessary to modify our approach and construct a basically new model as shown in 
Figure 21. This new model represents a change in mass attachment that complicates 
the coupling of mass to the column matrix of variables. 

Ki K2 Kg K4 
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The displacements of Mi.Mo and M 3 are the average of the two adjacent hinge 

, /yi ^ 2 \ 

displacements. Thus the dynamic force of may be represented by <jo^ I ~ — j . 


In order to expand the beam to a two dimensional model, the triangle concept is 
replaced with a massless rectangular pyramid as shown in Figure 22 and Figure 23, 





BOTTOM VIEW 


FIGURE 22. PYRAMID CONFIGURATION 


In the development of the matrix equation for this beam the mass displacement 
may be taken as average of two diagonal adjacent hinges; for convenience, the top left 
and bottom right of each pyramid base are taken. The displacement of spring is 


A = 
ai 




V + V > 

•'01 •^02 


and the displacement of 


A = 
a 2 



+ 




) 


The. first six lines of matrix equation (39) is the analytical model of the two- 
dimensional, four-spring beam. This model has no advantages for beam problems over 
equation 23 and is more complicated; but the corner hinges allow for modal coupling of 
a plate problem. 
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The plate made up of four-spring beams may be represented as in Figure 24. 



FIGURE 24. FOUR SPRING PLATE MODEL 

There are three horizontal beams and three vertical beams. The four corners of the 
plate are not a part of any beam. However, when the assumption that the pyramids are 
rigid is applied, the three hinges of the corner segments that are attached to beams will 
determine the displacement of the corners. 

Because of the coupling, the number of displacement variables of the plate is 
not equal to the sum of the six beams and the four corners. This fact reduces the 
number of required equations to make up a plate matrix equation. Each beam model 
must include all of the variables of that beam. When fewer equations are required, the 
fifth and sixth lines of the equation may be deleted and line seven substituted as the fifth 
equation of the beam set. The plate of Figure 24 will require a rectangular matrix of 
60 by 36 which reduces to 36 by 36 when the end conditions are inserted. Equation 40 
gives the analytical model of the plate. The end conditions are selected in the same 
manner as the single beams and all the variations as described for the beams may be 
applied to the plate matrix. 

The stiffness constant of the two dimensional beam model is 
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Equation (39) 
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The plate made up of four-spring beams may be represented as in Figure 24. 



FIGURE 24. FOUR SPRING PLATE MODEL 

There are three horizontal beams and three vertical beams. The four corners of the 
plate are not a part of any beam. However, when the assumption that the pyramids are 
rigid is applied, the three hinges of the corner segments that are attached to beams will 
determine the displacement of the corners. 

Because of the coupling, the number of displacement variables of the plate is 
not equal to the sum of the six beams and the four corners. This fact reduces the 
number of required equations to make up a plate matrix equation. Each beam model 
must include all of the variables of that beam. When fewer equations are required, the 
fifth and sixth lines of the equation may be deleted and line seven substituted as the fifth 
equation of the beam set. The plate of Figure 24 will require a rectangular matrix of 
60 by 36 which reduces to 36 by 36 when the end conditions are inserted. Equation 40 
gives the analytical model of the plate. The end conditions are selected in the same 
maimer as the single beams and all the variations as described for the beams may be 
applied to the plate matrix. 

The stiffness constant of the two dimensional beam model is 
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where v is Poissoai's Ratio and the other nomenclature is the same as in the beam. 

The 4 by 4 spring plate solutions will be accurate for only the first mode of the 
plate. Practical considerations would dictate a matrix up to the capacity of the available 
computer. For example a plate with 13 by 13 springs and good for 3 modes in each 
direction would require a 169 by 169 matrix; a plate good for 4 modes in each direction 
would require a 289 by 289 matrix. 

Since maximum amplitudes will be encountered in the lower modes, there is 
practical value from the use of the model. The Eigenvalue problems need only to be 
computed for a limited number of the lowest frequencies. Higher frequency solutions 
would have little meaning. 


SECTION VI. VARIATIONS OF THE PLATE MODEL 


There are two major variations of the plate model that have significance. There 
are adaptations of the model to allow the inclusion of stiffening beams, and the simula- 
tion of a cut-out plate for a complex beam structure. Both represent a use of additions 
or subtractions of springs and masses to the plate model. The matrix equation can be 
modified to meet the specific applications. 

Shell Model 


The analytical model can be modified to simulate a thin cylindrical shell struc- 
ture. The deflection is assumed in a radial direction. With r (p substituted into the 
matrix equation for i of the wraparound beam, the wraparound beams are continuous 
and therefore have no end conditions. This changes the number of unknowns and also the 
character of available equations to set up a suitable matrix. 

The adaptation of analytical principles used in this paper will be extended to 
cylinders and published at a later date. 




SECTION vn. CONCLUSIONS 


The dynamics of a beam or plate can be simulated by a model made up of mass- 
less constraints, springs, and concentrated masses. By unique arrangement of the 
mechanical elements, a mathematical model may be produced for structures not gener- 
ally solvable by usual techniques. The analytical approach has been derived and the re- 
sults proven by the use of simple hand-solvable models. Extension to more complicated 
structures is possible using digital computers for matrix Eigenvalue and Eigenvector 
solutions. 

George C. Marshall Space Flight Center 
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Huntsville, Alabama, February 1965 
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nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA aaivities but not necessarily reporting the results of individual 
NASA-programmed scientific eflforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 


SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washinston^ D.C. 20546 



